Manifestation of the spin-Hall effect through transport measurements in the 

mesoscopic regime 
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We study theoretically the manifestation of the spin-Hall effect in a two-dimensional electronic 
system with Rashba spin-orbit coupling via dc-transport measurements in realistic mesoscopic H- 
shape structures. The Landauer-Buttiker formalism is used to model samples with mobilities and 
Rashba coupling strengths of current experiments and to demonstrate the appearance of a measur- 
able Rashba-coupling dependent voltage. This type of measurement requires only metal contacts, 
i.e., no magnetic elements are present. We also confirm the robustness of the intrinsic spin-Hall effect 
against disorder in the mesoscopic metallic regime in agreement with results of exact diagonalization 
studies in the bulk. 



Introduction. The ability to manipulate electronically 
spins and to generate spin currents in semiconductors is 
the sine qua non for the full development of semicon- 
ductor based spintronics 0]. The control of spin and 
spin-currents without applying external magnetic fields 
can be achieved through the spin-orbit (SO) coupling, 
which acts as an effective momentum-dependent Zeeman 
field. Within this context, the newly discovered intrin- 
sic spin-Hall effect (SHE) in p-doped semicondcutors by 
Murakami et al. '5] and in a two-dimensional electron sys- 
tem (2DES) by Sinova et al. 3] offers new possibilites for 
spin current manipulation and generation in high mobil- 
ity paramagnetic semiconductor systems. The intrinsic 
spin-Hall effect represents a spin-current response gener- 
ated perpendicular to the driving electric field. The spins 
are tilted out of the plane due to the torque imparted by 
the SO coupling induced effective Zeeman field. In the 
Rashba SO coupled 2DESs the bulk intrinsic spin-Hall 
conductivity was found to have a value of e/87r in the 
clean limit for the case of both spin-split subbands being 
occupied and decreases linearly with the electron density 
for single spin-split subband occupation (3[. 

The SHE discovery has generated a tremendous inter- 
est in the research community |1I3,|E0,I1I1II3,IH 
in III III 111 El- Similarly to the long-standing debate 
on the origin of the anomalous Hall effect (AHE) [l^ . 
the robustness of the bulk intrinsic SHE against disorder 
and how it is related to the scattering mediated extrin- 
sic spin-Hall effect ^1 HolE ^ has been the focus of an 
intense theoretical debate 0lijl 111 III III il^] • While 
it was understood originally [3j that, unlike the quan- 
tum Hall effect, the universal value of the intrinsic SHE 
in the Rasbha SO coupled 2DESs will be reduced when- 
ever the disorder broadening is larger or similar to the 
SO coupling splitting, as was verified within a standard 
Born-approximation treatment jjj, taking into account 
the ladder vertex corrections through various methods 
suggests that the bulk spin-Hall conductivity vanishes in 



the weak disorder dc- limit H However, these re- 

sults have been challenged by other analytical calcula- 
tions which also consider ladder vertex corrections |l2l |. 
Other recent studies ,13], using arguments that echo the 
long-standing debate between skew and side-jump scat- 
tering in the AHE, have argued that the intrinsic SHE 
vanishes in all regimes. 

Given the ferraginous collection of analytical results, 
unbiased numerical calculations are needed to shed light 
on this controversy. An exact diagonalization treatment 
of disorder has shown that the bulk intrinsic SHE 
is robust against weak disorder. In addition, several nu- 
merical studies utilizing the Landauer-Buttiker (LB) for- 
malism in a tight-binding model representation of the 
Rashba SO coupling Hamiltonian in the presence of dis- 
order show similar conclusion in the limit which corre- 
sponds to the continuum effective mass model (see below) 

111 El H 113. 

In this paper we address a key question that has yet 
to be addressed directly: How to measure the intrin- 
sic spin-Hall effect through transport measurements? All 
previous numerical studies have focused on the contro- 
versy regarding the robustness of the effect against dis- 
order and how disorder, Rashba coupling strength, etc., 
change the continuum effective mass model value of c/Stt 
in the various models. Recently a new theoretical ques- 
tion has arisen, whether the dissipationless currents or 
spin background currents can lead to spin accumulation 
or to a steady signal which can manifest such an effect 
H . This question can be addressed unambiguously in 
the mesoscopic regime where the effect of the leads and 
disorder can be taken into account through the explicit 
treatment of voltage and current probes within the LB 
formalism 0, lisl \Wi |2]| . Here we consider an H-shape 
structure shown in Fig. ^ to demonstrate the appear- 
ance of the spin-Hall effect through dc-transport mea- 
surements without any magnetic elements. A current is 
driven from lead 2 to 1 in the lower leg shown in Fig.n(a) 
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and the rest of the contacts are voltage probes, i.e., leads 
with total zero current. Then, for typical current utilized 
in experiments and typical parameters we find a voltage 
difference AV34 = V3 — V4 dependent on the Rashba cou- 
pling and coupled to the spin-Hall conductance defined 
at, e.g., probe 6. The variation of this voltage, relative 
to the residual voltage obtained at zero Rashba coupling, 
is of the order iiV for the maximum system sizes that we 
can model (~ 0.1/im). 

\\aM 




FIG. 1: (a) Mesoscopic H-bar probe with metallic leads, (b) 
The continuum model is converted into a tight-binding model 
and the effects of the leads are treated exactly through a self- 
energy. 

In the mesoscopic regime, several of the controversies 
that have arisen from the study of the bulk spin-transport 
coefficients can be addressed. Within this regime the only 
assumption made in describing the transport through the 
sample via the current and voltage probes, other than the 
applicability of the tight-binding approximation for the 
SO coupled electronic structure, is that such contacts are 
perfectly metallic, i.e., an exact analytical expression is 
known for their Green's functions. 

Several numerical studies have addressed the robust- 
ness of the intrinsic spin-Hall effect within the mesoscopic 
regime utilizing the LB formalism in a tight-binding 
model representation of the Rashba SO coupled Hamil- 
tonian El E H EH . Xie et al. ^ obtained the ex- 
pected universal spin-Hall conductance in the weak dis- 
order limit with the Fermi energy, Ep, at the band cen- 
ter Ep = 0, but observed that spin-Hall conductance 
decreased rapidly with system size due to localization ef- 
fects. These results are in direct contradiction to the 
general notion that the SO coupled 2DESs exhibit a de- 
localized region [iH |2^ and, perhaps more importantly, 
that at the band center of this model anvHall coefhcient 
vanishes due to electron-hole symmetry jlj, |ig . Perhaps 
the most compendious of these tight-binding model nu- 
merical studies are Refs. Iisl and ITil where the expected 



symmetry of the Hall conductance with respect to Ep 
and the expected metal-insulator transition as a function 
of SO coupling and disorder strength . 

Model Hamiltonian and LB treatment of the spin-Hall 
effect. The experimental detection of the SHE through 
electrical means is conceptually challenging. Given the 
controversy surrounding the nature of the spin-currents 
generated by electric fields, a measurement of the voltage 
between two metallic contacts 01 appears to be the most 
promising dc-transport approach to unambiguously de- 
termine the presence of the SHE signal. We focus our 
attention then in the proposed H-shape device shown 
in Fig. n s-iid demonstrate that within the mesoscopic 
metallic regime the intrinsic SHE is exhibited through 
the change in a voltage difference between two contacts 
as the Rashba spin-orbit coupling is varied. 

The continuum effective mass model described by the 
2DES Hamiltonian with the Rashba SO interactions is 
given hy H = + X{axPy - ^yPx) + ffrf,:,,, where the 
second term is the Rashba SO coupling [23, '^^^ ^o 
the asymmetry of the confining potential and H^ig is 
the disorder potential To model the complex geometry 
of our disordered conductor within the LB formalism we 
use the tight-binding (or finite differences) approxima- 
tion [ij, |2lj . Within this approximation the continuum 
effective mass envelope function Hamiltonian becomes: 

j 

+ Xl^slt^J+'^-i ^ Su'^j+'^-t) + ^-c-], (1) 

i 

where t — /2m*af^ and tso = -^/2ao, oq is the mesh 
lattice spacing, and 5 — ±aoi, ±aoy. The first term rep- 
resents a quenched disorder potential and disorder is in- 
troduced by randomly selecting the on-site energy Sj in 
the range [-W/2,W/2]. Within the leads the SO cou- 
pling is zero and therefore each lead should be consid- 
ered as having two independent spin-channels. These 
leads constitute reservoirs of electrons at chemical po- 
tential ^1, . . . , ^jv where N is the number of leads which 
we consider to be either four (lead 1-4 in Fig. ^ or six 
(leads 1-6 in Fig.^. 

In the low temperature limit ksT << Ep and for 
low bias-voltage, the particle current going through a 
particular channel is given within the LB formalism by 
tm Ip,a = {e/h)J2ga'Tp.a-q,a'[Vp " Vg], whcrc p labels 
the lead and Tp^a-,q.<j' is the transmission coefficient at 
the Fermi energy Ep between the (p, a) channel and 
the (g, a') channel. This transmission coefhcient is ob- 
tained by Tp^cr;q,a-' = ^^[^p,aG'^rq^^i G^] whcrC Tp^a- IS 

given by Tp^^iiJ) = i[i:^,^{i,j) - S^,^(z, j)], and the re- 
tarded and advance Green's function of the sample G^/"* 
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with the leads taken into account through the self energy 
T.plJ^{i,j) is given by 

Here the position representation of the matrices Tp,a, 
G^, Hi J, and are in the subspace of the sample, 
i.e., it only includes sites in the sample. Since the 
SO coupled Hamiltonian preserves time reversal sym- 
metry, the transmission coefficients obey the relation 




like (In,Ga)As. We take the Rashba parameter A in the 
range from to 80 meV nm which are easily obtained in 
experiments [i^, El] , and we choose the electron concen- 
tration of approximately n2D = 10^^ cm~^. The Fermi 
energy is obtained from the chosen electron concentration 
assuming an infinite 2DES. This gives a small difference 
of a few percent when considering our finite tight-binding 
model but the leads themselves are the ones providing 
the reservoir of electrons and therefore such fluctuations 
are small as verified by direct numerical calculations (not 
shown) . 

In Fig. 2 we show the spin Hall conductance Gsh as a 
function of Rashba parameter for the H-shape structure 
when current flows through the bottom leg as indicated 
in Fig. n Here we consider the system with 6 leads and 
with leg lengths L varying from 90 to 140 nm. The total 
size of the system is L in the a:-direction and L/2 in 
the y-direction. The horizontal connection bar is L/6 by 
L/6. The width of the legs is L/6 with the attached leads 
of the same width. These ratios were chosen for typical 
fabricated samples (of larger system size) but any shape 
is feasible to do and a search for an optimal geometry is 
underway |29| . 
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FIG. 2: Spin-Hall conductance defined at lead 6 (shown in 
Fig. vs. spin-orbit coupling strength for different size sys- 
tems for mo = 0.05me,/i = 250, OOOcm^/V^s and flowing cur- 
rent of lOnA in the bottom leg. Here L is divided in 42 points. 
Only a few disorder realizations are needed for convergence 
in samples with these mobilities. 

Within the above formalism the spin current through 
each channel is given by ,^ = (e/47r) Y,qo' Tp,a-q,a' [Vp - 
Vq] and through this we define a spin-Hall conductance 
as 

Gsh = (3) 

Vl — V2 

as indicated in Fig. ^ All the voltages are obtained by 
imposing the boundary conditions + =0 for i=3 
through 6, + /i| = 1 and /2.1 + /2X = —1- The 
arbitrary zero of voltage is fixed by setting V2 = 0. These 
are later translated to a realistic voltages by setting the 
current to a typical value of 10 nA. 

Results and discussion. In order to address the key is- 
sue of how the spin-Hall effect can manifest itself through 
a dc-transport measurement without ferromagnetic con- 
tacts we choose realistic parameters for our calculations 
which model currently available systems [2^. We con- 
sider an effective mass of m* — O.OSwe and a disorder 
strength oiW = 0.09 meV corresponding to the mobility 
of 250,000 cm^/Vs, which is typical for a semiconductor 
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FIG. 3: The voltage difference between leads 3 and 4 as a 
function of Rashba couphng for H-probe for different size sys- 
tems and meshes for mo = 0.05me, fj, = 250, OOOcm^ /Vs and 
flowing current of lOnA in the bottom leg. 

This particular H-shape structure allows for minimal- 
ization of the residual voltage drop due to charge flow 
which is of the order of few hundreds of nanovolts. Gsh 
is calculated accordingly to Eq. |31 The calculations are 
conducted for a few different meshes Ni — L / qq to check 
the convergence of results for oq ^ 0. The magnitude 
of Gsh is around 0.2-0.6 in e/Sir units for Rashba cou- 
pling 70-80 meV nm and L = 120 - 140 nm. We also 
note that within these parameters we are well within the 
metallic regime |lH l22j and both spin-split subbands are 
occupied. 
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The spin Hall conductance cannot be measured by the 
paramagnetic leads and ferromagnetic leads can intro- 
duce spurious effects coming from the impedance mist- 
match preventing ballistic contacts |25ll3Ctl3l| . However, 
we expect that the spin current which flows between leads 
5 and 6 can generate a secondary effect, the induction of 
a voltage difference in the top leg between leads 3 and 4. 
This is in the same spirit as the initially proposed set-up 
by Hirsch but in a far simpler configuration with- 
out the need of considering a bridged conductor nor any 
unknown scattering mechanisms other than the effects 
of disorder which is actually small in this case [isl . 
We illustrate this in Fig. 3 where we show the nonzero 
voltage difference AV34(A) — AVi34(A = 0) as a function 
of Rashba coupling for different size systems, different 
meshes and with and without the additional leads 5 and 
6. We find the increase of AV34 with the increase of 
Rashba coupling. The induced voltage variation is of the 
order of few /iV for A — 60 — 80 meV nm and can be easily 
measured. We also note that the inclusion or omissions of 
leads 5 and 6 which in reality cannot measure directly the 
spin Hall conductance and do not influence the voltage 
difference; they do however influence the total residual 
voltage background AV34(A = 0). The convergence of 
the results as a function of mesh size is illustrated in Fig. 
3 which show the results for the larger meshes calculated. 

We also note that at such disorder strength the results 
do not depend on the disorder strength if we, say, double 
it. This indicates two important facts: one, that the 
voltage induced is originating from an intrinsic effect, 
and two, that the effect is not related to the constant in- 
plane polarization induced by a current in a 2DES with 
Rashba coupling, the Levitov effect ji^l , since this effect 
is proportional to the mobility. 

Summary. We have calculated, as a function of the 
Rashba SO coupling strength A, the voltage drop AV34 
that occurs in an H-shape sample in a response to a driv- 
ing dc-current between leads 1 and 2. The voltage dif- 
ference closely follows the changes of the calculated SHE 
conductance, Gsh, with A. Moreover, AV34 increases 
with the increase of Rashba coupling for constant disor- 
der strength W which is a clear evidence that the ob- 
served voltage signal is directly connected with the in- 
trinsic SHE. Our work provides another confirmation of 
the robustness of the intrinsic SHE against weak disor- 
der and shows the feasibility of detecting SHE signals 
through dc-transport measurements in structures with 
realistic experimental parameters. 
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